function Test()

%[p, vt, Phi, Delta, S] = PDEFD(0.2, 0.1, 50, 50, 0.005, 0.001, 1, 'CNS', 'put', 500, 40); 
%for i = 1 : length(S)
%	B(i) = BlackScholes(S(i), 0.2, 1, 50, 0.1, 'put');
%end
%
%[p, vt2, Phi, Delta, S2] = PDEFD(0.2, 0.1, 50, 50, 0.005, 0.001, 1, 'CNS', 'put', 500, 15); 
%for i = 1 : length(S2)
%	B2(i) = BlackScholes(S2(i), 0.2, 1, 50, 0.1, 'put');
%end
%
%[p, vt3, Phi, Delta, S3] = PDEFD(0.2, 0.1, 50, 50, 0.005, 0.001, 1, 'CNS', 'put', 100, 40); 
%for i = 1 : length(S3)
%	B3(i) = BlackScholes(S3(i), 0.2, 1, 50, 0.1, 'put');
%end
%
%[p, vt4, Phi, Delta, S4] = PDEFD(0.2, 0.1, 50, 50, 0.005, 0.001, 1, 'CNS', 'put', 50, 10); 
%for i = 1 : length(S4)
%	B4(i) = BlackScholes(S4(i), 0.2, 1, 50, 0.1, 'put');
%end
%
%save fdcnsput.mat -V6 
%
[p, vt, Phi, Delta, S] = PDEFD(0.2, 0.1, 50, 50, 0.005, 0.001, 1, 'CNS', 'call', 40, 60); 
for i = 1 : length(S)
	B(i) = BlackScholes(S(i), 0.2, 1, 50, 0.1, 'call');
end

[p, vt2, Phi, Delta, S2] = PDEFD(0.2, 0.1, 50, 50, 0.005, 0.001, 1, 'CNS', 'call', 30, 70); 
for i = 1 : length(S2)
	B2(i) = BlackScholes(S2(i), 0.2, 1, 50, 0.1, 'call');
end

[p, vt3, Phi, Delta, S3] = PDEFD(0.2, 0.1, 50, 50, 0.005, 0.001, 1, 'CNS', 'call', 30, 90); 
for i = 1 : length(S3)
	B3(i) = BlackScholes(S3(i), 0.2, 1, 50, 0.1, 'call');
end

[p, vt4, Phi, Delta, S4] = PDEFD(0.2, 0.1, 50, 50, 0.005, 0.001, 1, 'CNS', 'call', 40, 50); 
for i = 1 : length(S4)
	B4(i) = BlackScholes(S4(i), 0.2, 1, 50, 0.1, 'call');
end


save fdcnscall.mat -V6 
